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Abstract 

We study the motion of a classical point body of mass M, mov- 
ing under the action of a constant force of intensity E and immersed 
in a Vlasov fluid of free particles, interacting with the body via a 
bounded short range potential ^f. We prove that if its initial velocity 
is large enough then the body escapes to infinity increasing its speed 
without any bound (runaway effect). Moreover, the body asymptoti- 
cally reaches a uniformly accelerated motion with acceleration E/M. 
We then discuss at a heuristic level the case in which \&(r) diverges 
at short distances like gr~ a , g,a > 0, by showing that the runaway 
effect still occurs if a < 2. 
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1 Introduction 



In the present paper we study the longtime behavior of a classical point body 
moving under the action of a constant force E and immersed in an infinitely 
extended medium. This problem has largely studied in the framework of 
kinetic theories, while here we want to obtain some rigorous results in the 
case of a fully Hamiltonian system. 

Of course, when the medium is absent the body performs an uniformly 
accelerated motion. In the presence of a medium the motion is much more 
complicated and its asymptotic behavior depends on the nature of the back- 
ground. In general, when the medium is homogeneous two different effects 
can happen: for strong body/medium interactions the motion of the body 
converges (in average) to an uniform one, while for weak body/medium in- 
teractions the body with a large initial velocity escapes to infinity increasing 
its velocity without bound. In the first case we have a reasonable model of 
viscous friction, in the second case we are in presence of a runaway particle 
(see [27] for a discussion of this effect in the framework of kinetic theories). 

Concerning the assumptions on the medium, a natural choice is to con- 
sider a system of infinitely many particles, mutually interacting via a pair 
potential $, and interacting with the body via a potential Unfortunately, 
this problem is very hard. Actually, the same definition of the model is not 
obvious, due to the difficulty to prove the existence of dynamics for infinitely 
extended systems. We recall that the time evolution of systems with infinitely 
many particles has been investigated in several papers: we only quote some 
of the main results [3,11,17,19-21,24,25], other references can be found 
therein. Long time estimates drastically depend on the spatial dimensions of 
the system, and good enough estimates are known only for one dimensional 
models with bounded interactions [10]. The runaway effect has been shown 
to occur for large external force E when the medium is confined in a tube [6], 
and for any strength of E in the case of particles moving along a line [7]. For 
singular interactions, the strict one-dimensional case becomes too particular, 
while in general we can only make conjectures. In [6], for interaction singular 
as r~ a the threshold case is conjectured to be a = 2, but a rigorous proof 
seems too hard. 

To make further steps we must simplify the medium. A possible way 
is to consider a mean field approximation, i.e. to put the body in a Vlasov 
fluid. In the mean field approximation, well known in Astronomy and in 
Plasma physics, the background is composed by a gas of particles of mass 
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going to zero and the number of particles per unit volume going to infinity 
in such a way that the mass density remains finite. The mean field limit 
of the dynamics of interacting particle systems in the case of finite total 
mass has been studied in [4,18,28,30] for bounded interaction and in [22] for 
singular force like r~ a , a < 1. The limit in the case of infinite total mass has 
been proved in [5] in the case of bounded interaction and spatial dimension 
one. For unbounded mass, the existence of the dynamics has been proved in 
three dimensions of a Vlasov fluid with bounded interactions [9], and in two 
dimension for the Helmholtz equation [12]. 

In this framework some results have been recently obtained. In [13, 15] 
the interaction \& is assumed of hard core form, precisely ^(Irl) = oo for 
r e A, ^ = otherwise, where A is a convex set of 1R 3 (this means that 
an element of the medium evolves freely out of A and interacts elastically 
on <9A). The initial phase space density of the medium is assumed spatially 
homogeneous out of A and with Gaussian distribution of the velocities. It 
is proved the existence of a stationary state in the motion of the body for 
any intensity of a spatially constant force E, and a detailed analysis of the 
approach to the limiting velocity is given. This result has been extended to 
forces not constant in space in [14], and a similar result has been obtained 
changing the elastic boundary condition with the diffusive ones (see [1,2] for 
analytical and numerical results). In the case of singular interaction of the 
form r~ a , a > the problem is open. We only recall the paper [8], where the 
runaway particle effect has been proved when \1/ has the form r~ a , a < 2, the 
medium is composed by a Vlasov fluid, the mutual interaction $ is bounded, 
the influence of the body on the fluid motion is neglected, and the system 
has initially an one-dimensional symmetry. 

In the present paper we assume $ = 0, the body/medium interaction ^ 
bounded and with a short range, the background initially at rest, and we 
find not only the runaway effect but also the exact asymptotic velocity of the 
body. Then we discuss at an heuristic level the case in which \I/ is always at 
short range but diverging at short distances as r~ a , a < 2, case in which also 
the runaway effect is expected to take place. In [6] it has been conjectured 
that a = 2 is the threshold value for the existence or not of the runaway 
effect. A discussion on the existence of a stationary state for this model [16] 
improves this conjecture. 

The plan of the paper is the following. In Section 2 we introduce the 
model and state the main result. In Section 3 we prove the result for \l/ 
bounded. Finally, in Section 4 we discuss at a heuristic level the case of 
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singular interaction. 



2 The model 

We consider the motion of a point body of mass M and position £, under 
the action of a constant force E of intensity E and directed along the x-axis, 
i.e. E = (E,0,0). The body is immersed in a fluid of free particles, that 
interact with the body via a force of pair potential \l/(|r|). We assume the 
medium be a Vlasov fluid of free particles in three dimensions, namely the 
pair of functions, 

(x, v) -> {X(x, v; t), V(x, v; t)), f Q (x, v) -> f(x, «; £), (2.1) 

with (cc, t>, f) 6 I 3 x I 3 x K, solution to 

X(aj,»;t) = V(x,«;t), (2.2) 

V(aj,i>;*) = -V*(|X(oj, -^(t)|), (2.3) 

(X(aj,t;;0),V(aj,»;0)) = (a:,t»), (2.4) 

/(X(x, w; t), V(x, v; t),t) = f (x, v), (2.5) 

where £(t) is the coordinate of the point body, / is the mass density in phase 
space of the free gas, and (X, V) are the characteristics of the fluid. 

In the case of smooth initial data, f(X(x, v; t), V(x, v; t),t) satisfies the 
differential equation, 

(d t + v-V x - W(|oj - £(£)|) ■ V v )f(x, v- 1) = 0. (2.6) 

The system f l2.2|) - fl2.5p is completed by the equation of motion of the body, 

M'i(t) = E- JdxdvVt*{\t{t)-x\))f(x,v;t), (2.7) 
with boundary condition 

m=to, m = t - (2.8) 

This kind of system, in which a Vlasov fluid is coupled to a massive body, 
has been firstly introduced in connection with the so-called piston problem, 
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see [23] and also [26] and references quoted therein. In our problem the fluid is 
unbounded and so, in principle, the existence of the solutions is not obvious. 
In fact, it is easy to exhibit initial conditions for which the time evolution 
produces singularity in the motion after a finite time. This happens because 
very far away particles could arrive very fastly close to the body. These 
situations are pathologic from a physical point of view and can be removed 
in the case of bounded interaction by assuming that initially 

f ( x , v )< p (J-y 2 e-^ 2 /\ p >0, (2.9) 

where (3 = (kT)^ 1 > 0, T is the temperature, and k the Boltzmann constant, 
i.e. / is bounded by a homogeneous Gibbs state. In this case the proof of 
the existence of the solutions is quite easy. 



3 Bounded interaction 

Assume ^(H) be a twice differentiate function of r 6 R 3 and let exist a 
positive constant r < oo such that \I/(r) = if r > r . We assume that 
initially the point body is in the origin, £ = (0,0,0), with velocity along 
the x-axis, £ = (£o,0, 0), £o > 0, and the fluid is at rest with constant 
density p(x, 0) = po > 0. By symmetry the body will move along the x-axis, 
i.e. £(t) = (£(£), 0,0) for some f(t) G R. 

Theorem 3.1 For each intensity E of the force, there exists a threshold 
valued such that for £o > £,* the body escapes to infinity with (asymptotically) 
a uniformly accelerated motion, 

Urn = * (3.1) 

Proof. We sketch the proof by emphasizing the main steps. We need to prove 
that the friction force exerted on the body by the fluid is initially bounded 
and vanishes asymptotically in time. For this purpose we study the time 
evolution of a particle (i.e. a characteristic) of the fluid. 

First, we remark that if is large enough then each fluid particle in- 
teracts with the body at most one time. To prove this property we argue 
in the (non-inertial) reference frame where the body is at rest. During the 
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collision, the fluid particle of coordinates x(t) = x 2 (t), x 3 (t)) feels two 

forces, — V\l/ and —£(t) = (—£(£), 0,0). Assume for the moment that 



and 



sup |£(t)| < A < oo, (3.2) 

t>o 



£(*)>§ Vt>0. (3.3) 



(In the sequel, these properties will be proved to be verified for A and £o 
large enough). 

Denote by r the time at which the interaction begins and by r + 8 the 
time at which it finishes. At time r the velocity of the fluid particle is 
x(t) = — £(t) and during the interaction its absolute value decreases at 
most by (B + A) 8, where B = sup f , |\I/'(r)|. Hence 



that for large £ implies 



8 < -. ^ , (3.4) 

-\t{r)\-25(A + B)' 



8 < S^. (3.5) 



Hence, for large £o the motion of the fluid particle does not change its direc- 
tion and the particle escapes to — oo. 

This discussion suggests also the proof of the assumption (13.21) : the body 
could interact only with the fluid particles that are present in a cylinder of 
finite length and so the acceleration of the body is bounded uniformly in 
time. Hence (13.21) holds for A large enough. 

Now, we come back to the reference frame in which the medium is initially 
at rest and we evaluate the momentum variation of the fluid particle during 
the collision with the body. Actually, a bound on \x(r + 8)\ can be easily 
achieved by (13. 5p . getting |jb(r + 5)| < const. H^'Hoo |£( r )|~\ but we look for 
a more sophisticated bound useful in the next section. In fact, we expect that 
the momentum variation is really very small for large velocities: the effects 
produced when x\ > £ compensates in part those produced when x\ < £. To 
show this fact we introduce the quantity p(t) = (pi(t) , p 2 (t) , p 3 (t)) defined as 

« t)=m + *MLz^ m - im (3.6) 
\x(t)-£(t)\ 2 
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Obviously, before and after the collision, 

pit) = x(t). 



(3.7) 



We shall show that during the collision p(t) is quite small and consequently 
the momentum variation is small. Recall [r, r + 5} denotes the time interval 
during which the collision occurs. By (13. 7p and remembering that before the 
collision the fluid particle is at rest, we have, 



which implies 



x(t + 5) = 

\xi(t + 8)\ < 
\x 2 (t + 5)\ < 
Mr + 5)\ < 



ds p(s), 



ds \pi{s) 



-+5 



ds \p 2 (s) 



T+S 



ds |p 3 (s) 



(3.8) 

(3.9) 
(3.10) 
(3.11) 



We next use these equations to obtain finer estimates on \x(t + 6)\. From 
now on C denotes a generic positive constant whose numerical value may 
change from line to line. We start by studying \x 2 (t + S)\. From (13.60 we 
have, 



x 2 



1 M\x-£\) 



\x — £|- 

2x 2 *(\x-S\)(*-£)-(*-Z) 
\x-^\ 4 



+ x 2 ; ~7~r. ~ + 



x — £| \x — £| 



(3.12) 



We observe that \x 2 \ < \^'{\x — £|)| (by the equations of motion) and, for £ 
large enough, £ > ~|± — £|. By (13.101) . 

\x 2 (t + 5)\ < C j^ds \*X\x{s)-t(s)\)\ [1+ I^IU^s)- 2 ]. (3.13) 
Changing variable ds = £ _1 d£ into the last integral and observing that 

^S V(\x(s)-£( S )\)<C\\nocC\ 
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where 



we have 



C = min J(s) \ 



\x 2 (r + S)\ < C 



\^\ 



c 



1^1 



C 2 



(3.14) 



(3.15) 



Due the cylindrical symmetry the same estimate holds for \x 3 (t + 5)\. We 
then turn to the estimate of |xi(r + <5) | . By (I3.6P and the equations of motion, 



n\*-t\) 

\x — £\\x — £| 



*(!*-*!) 

\x-k\ 2 



(xi - g)(x\ + ±1) + {£- x 1 )(x 2 x 2 + x 3 x 3 ) 
2(i-x 1 )(£-x)-{£-x 



(3.16) 



By plugging the bound (13 . 15[) and again equations of motion into (|3.16p . by 
and recalling the definition (13.141) we obtain 



|ii(r + 5)|<cMk 



1 + 



1^1 



C 2 



+ 



c 4 



+ 



(3.17) 



Hence, since H^Hoo < °°; each fluid particle increases the ^-component of its 
momentum by a quantity at most proportional to (~ 3 . Since the body meets 
per unit time a number of fluid particles proportional to £(£), it looses per 
unit time an amount of momentum bounded as below, 



[lost momentum per unit time] < C £(t) 2 . 



(3.18) 



The other steps in the proof are easy: we show that, for Ci,£o > large 
enough, 



Mi(t) > ^io + Et- 



o 



(3.19) 



Indeed, denote by t* the supremum of the times T such that the inequality 
(ETUI is verified for any t E [0,T]. Since Mf(0) = M£ , we have t* > by 
continuity. We decompose the time interval [0, £*) into a partition of intervals 



5k of length of the order <p y Yli=i ^* J j where 4>{i) is the solution to the integral 
equation 

Ci 



M(f>(t) 



3M • 

-r- £o + Et- 



ds 



Since ( 13.191) suddenly implies the assumption (13. 3p for £ large enough, the 
bound (I3.18P is valid in each interval 8^. It follows that the approach to 
the asymptotic (uniformly accelerated) motion of £(t) is faster than that of 
4>{t), whence t* = +oo for large enough C\ and £o- So (I3.19P is valid at any 
positive time and the theorem is proved. □ 

Remark 3.2 We have assumed the density of the fluid constant in the whole 
space. Of course, the result and the proof do not change if we assume the 
initial density po(x) with an axial symmetry around the x-axis and po{x) — > 
const, fast enough as \x\ — > oo. Indeed, this only changes the transient 
evolution, but the asymptotic behavior of the motion remains unchanged. 

It is natural to ask what happens if we assume the fluid initially dis- 
tributed as a Gibbs state like (I2.9P (Theorem 13.11 here corresponds to the 
case P = oo). For E larger than the result is still true [29], while for 
smaller E it depends on the sign of \l/'(r): we conjecture that if \?'(r) > 
the result remains valid; while, if \&'(r) < for some r, the quantity £(t)/t 
should converge to some constant possibly smaller than E/M. In fact, the 
body could capture a small part of fluid, which remains with it forever, with 
the only effect to increase the effective inertia (i.e. mass) of the body. 



4 Singular interaction 

The spirit of this section is different from the previous one, where we have 
given a mathematically rigorous proof on the asymptotic behavior of the 
system. Here we confine ourself to a heuristic discussion of the main issues 
in a possible proof. 

We assume \I/(r) be a twice differentiable function for r > and there 
exist two positive constants ri,r < oo such that ^ = gr~ a ,g > 0, a > if 
r < r±, *Sf equals to a not increasing function if r\ < r < vq, \& = if r > tq. 
Initially we put the point body in the origin with a fix velocity £(0) = £o > 0, 
and the fluid at rest with a constant density p(x, 0) = p > 0. We discuss 
the following proposition, analogous to Theorem 13. 1[ 

Proposition 4.1 For each intensity E of the force, there exists a threshold 
valued such that for £ > the body escapes to infinity with (asymptotically) 
a uniformly accelerated motion, 

lim M = £. (4.1) 
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The strategy of the proof is similar to that used in the bounded case: we 
assume a behavior of the body and we evaluate the viscous friction produced 
by the fluid. Using an adiabatic invariant we show that the friction is very 
small and vanishes at long times, so that the behavior of the body is better 
than the assumption. Unfortunately, we do not prove some (minor) steps 
and the discussion remains only heuristic. 

In the present analysis two new points arises: a) the interaction is singular 
and so the particles of fluid near the x-axis drastically change their velocity 
whatever large is the velocity of the body; b) some particles of the fluid 
can collide infinitely many times. The first difficulty is solved by proving 
the these bad particles are very few when the velocity of the body becomes 
large; the second difficulty is solved by showing the the effect of recollisions 
is negligible. 

As in the bounded case we assume that 

sup \({t)\ < A < oo. (4.2) 

t>o 

We do not prove it rigorously. However physically it is quite obvious; more- 
over, at the end of the discussion it will be clear that a proof using this 
assumption should imply that the acceleration goes to zero in average on a 
time interval of length vanishing as t — > oo. But a rigorous proof of assump- 
tion (14. 2 p appears cumbersome and not trivial. 

We explicitly discuss the more difficult case a > 1, starting with the 
study of the first collision. 

We analyze the motion of a fluid particle. We denote by rj its impact 
parameter, that is the distance between the incoming fluid particle and the x- 
axis. In the plane (22, £3) we consider the disks centered in the origin and 
of radius r\ k = 2 k r] , k = 0, 1, 2, . . . We choose r/ = £(t)~( 1+e \ e > 0, so that 
the particles with 77 < i]q contribute at most as C£(t) _2e in the momentum 
transferred to the body (per unit time). In fact, the area of this disk is 
7t?7q, the body meets per unit time a number of fluid particles proportional 
to £(£), in each collision it looses a momentum 2£(£) (corresponding to an 
elastic collision with a disk) and so the lost momentum due to particle hitting 
in this circle is at most C^ 2 {t)r]Q = C£(t)~ 2e . 

We now evaluate the effect of the collisions in the /c-annulus Dk \ -Dfc-i- 
First we remark that, by arguing as at the beginning of the previous section, 
it is possible to show that the body has only one collision with each fluid 
particle posed in a fc-annulus, and then it overcomes such particles. Hence, 
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the lost momentum is given by (13.171) . where now H^Hoo — C {i]k)~ a and the 
dominant term is II^H^,^) -3 (here we approximate (, see (I3.14p . by £(£)). 
The area of the /c-annulus is Crjl and the intensity is of the order of 
£ (t) . In conclusion, the contribution / to the lost momentum due to the fluid 
particles posed in all the annuli is bounded as 

k* 

I < c^) [(2a - 2)(1+e) ~ 21 J22 2k (~ a+1 \ Vk * > r„. (4.3) 

k=l 

If we choose e < (2/a) — 1 we obtain 

I<Ci{ty 2 \ (4.4) 

that is a bound analogous to the one obtained for the first disk. 

We now roughly discuss the effect of multiple collisions. For the moment, 
we suppose that the body is governed by a uniformly accelerated motion. 
The only possibility for a fluid particle to hit the body many times is to be 
quite close to the x-axis. In particular, we start by evaluating how much 
close for hitting the body twice. Denote by a the impact parameter of the 
fluid particle and let t\ be the time at which the first collision occurs. Due 
to the monotonicity of the potential, during the collision the fluid particle 
feels a repulsive force from the which varies its orthogonal velocity by 

a quantity 

rti+S 

5v ± ~ ds\V$(x(s)-£(s))\—. (4.5) 
Jti r o 

Changing variable ds = £ -1 d£ we obtain 

Sv^ZC-Hhy'r-*, (4.6) 

where r m is the minimal distance between the fluid particle and the body. 
Since r m fa C £~ 2 / a we conclude that 

5v L >Cai(h). (4.7) 

After the collision, the fluid particle has gained a velocity which is double 
than the speed of the body. But the particle then preserves its velocity while 
the body accelerates almost uniformly, and so it can be reached again after 
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a time of the order of In conclusion, the fluid particle can hit the body 

twice if its impact parameter a is very small, 

a<Ce(ti)- 2 . (4.8) 

Actually, we are interested in the case when the fluid particle has a third 
collision. For this to happen the impact parameter has to be very small, 
precisely 

^<C^ir 2 ^ 2 )- 2 <C^i)- 4 (4.9) 

where t 2 is the time of the second collision and we used £(£2) > 

To apply these ideas to the true case, some modifications are needed: 
taking into account that the motion of the body is not exactly uniformly 
accelerated, but with a small perturbation; moreover, to evaluate the effect 
of multiple collisions at a fixed time, we should go back in the past up to the 
initial time; finally, for the repulsive force we need a lower bound (not an 
approximation as in previous discussion). We do not perform this analysis in 
detail, but we hope to have convinced the reader that the effects of multiple 
collisions are negligible. 

The further steps are similar to that of the bounded case: we assume that 

q /\/f ft n 

Mat)>—Zo + Et- / dsj— , (4.10) 

4 Jo £(s) e 

and by ( 14.4ft we prove that the approach is faster. Then the previous bound 
holds at any positive time and the body motion converge to a uniform accel- 
erated one as t — > 00. 
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